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ABSTRACT 

The  location  of  Iraq  within  the  mid  latitudes  made  it  one  of  the  places  which  receive  large  quantities  of  solar 
radiation.  Due  to  the  global  trending  of  investing  this  energy  and  for  environmental,  healthy  and  economical 
considerations,  there  are  great  opportunities  to  construct  projects  that  invest  the  solar  energy  to  generate  power  in  Iraq, 
which  requires  providing  models  specialized  in  locating  the  spatial  and  temporal  distribution  of  the  incident  solar  radiation 
on  the  inclined  surfaces  to  provide  more  solar  energy,  In  this  study  a mathematical  program  has  been  proposed  by  using 
MATLAP  by  using  two  kinds  of  models  isotropic  and  an  isotropic  to  locate  the  spatial  and  temporal  distribution  hourly 
solar  radiation  incident  on  inclined  surfaces  at  latitude  (25°-  40°)  and  longitude.  (35°-  50°),  considering  the  inclination 
angle  is  the  optimum  angle  at  all  times.  It  shows  the  possibility  of  having  more  solar  energy  on  the  inclined  surface 
compared  with  the  horizontal  surface;  the  difference  was  obvious  in  the  months  (January,  February,  March,  October, 
November  and  December).  The  models  mentioned  above  were  also  used  in  determining  the  spatial  and  temporal 
distribution  of  the  incident  solar  radiation  on  the  inclined  surface  under  different  sky  conditions  (clear,  semi-cloudy  and 
cloudy)  above  Iraq  and  the  neighboring  countries,  it  indicates  that  the  difference  in  the  solar  energy  between  the  three 
conditions  is  larger  in  January  and  December  compared  with  other  months. 
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INTRODUCTION 

In  developing  countries,  such  as  Iraq,  interest  in  solar  energy  applications  has  been  growing  in  providing 
electricity  and  water  supply  in  different  areas.  Understanding  solar  radiation  data  is  essential  for  modeling  solar  energy 
systems.  Solar  radiation  is  used  directly  to  produce  electricity  for  photo  voltaic  (PV)  systems  and  solar  thermal  systems. 
Therefore,  precise  knowledge  of  historical  global  solar  radiation  at  a location  of  study  is  required  for  the  design  and 
estimation  of  the  performance  of  any  solar  energy  system.  Knowledge  of  solar  radiation  incident  on  an  inclined  surface  is 
essential  to  architects  and  engineers  for  energy-efficient  building  designs  and  solar  energy  applications.  The  global 
radiation  incident  on  a tilted  surface  consists  of  components  direct,  diffuse  and  reflected  from  the  ground.  On  a hourly 
database,  the  direct  radiation  can  be  calculated  by  geometric  projections  (ratio  of  the  incidence  angle  to  the  solar  zenith 
angle).  The  reflected  radiation  has  a small  effect  on  calculations  and  may  be  calculated  with  an  isotropic  model.  Both 
components  present  dependence  of  measures  in  incidence  or  horizontal  surface  the  meteorological  stations  measure  global 
and  diffuse  solar  radiation  intensities  mostly  on  horizontal  surfaces  only  [1]  Whereas,  the  stationary  solar  conversion 
systems  (both  for  the  production  of  electrical  and  thermal  energy,  i.e.  photo  voltaic  cells  and  flat  plate  collectors)  are  tilted 
towards  the  sun  in  order  to  maximize  the  amount  of  solar  radiation  incident  on  the  collector  or  cell  surface.  But,  the 
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availability  of  these  data  on  tilted  surfaces  is  very  rare;  therefore,  the  tilted  surface  radiation  in  most  cases  is  calculated 
from  measured  global  horizontal  radiation  [2,  3]. 

One  type  of  model  estimates  beam  and  diffuse  components  from  global  horizontal  radiation  and  the  other  predicts 
tilted  surface  radiation  [3]  this  is  usually  an  inclined  surface,  in  order  to  maximize  the  solar  radiation  received,  whether  by 
using  a fixed  inclined  mounting  position  or  a sun  tracker  system.  [4]  The  total  radiation  on  a tilted  surface  consists  of  three 
components:  beam,  reflected  radiation  from  the  ground  and  diffuse  from  the  all  part  of  the  sky.  The  direct  and  reflected 
components  can  be  computed  with  good  accuracy  by  using  simple  algorithms,  but  the  nature  of  diffuse  part  is  more 
complicated.  Calculation  of  diffuse  radiation  requires  information  of  both  global  and  direct  radiation  incident  on  a 
horizontal  surface  at  the  same  time  [5]  Solar  radiation  on  the  earth  can  be  classified  as  Direct  radiation  this  radiation  comes 
directly  from  the  sun  without  any  change  in  its  direction,  Diffuse  radiation:  This  radiation  comes  from  all  over  the 
atmosphere  as  a result  of  reflection  and  scattering  by  clouds,  particles  in  the  atmosphere,  dust,  mountains,  trees,  buildings, 
the  ground  itself,  and  so  on,  Global  solar  radiation  consists  of  direct  and  diffuse  solar  radiation,  The  radiation  that  reflects 
from  surroundings  (so-called  albedo)  is  of  importance  for  some  surfaces  that  are  inclined  under  some  angle  to  the 
horizontal  surface  [7]. 

SOLAR  RADIATION  ON  INCLINED  SURFACE 

When  solar  radiation  enters  the  Earth’s  atmosphere,  a part  of  the  incident  energy  is  removed  by  scattering  or 
absorption  by  air  molecules,  clouds  and  particulate  matter  usually  referred  to  as  aerosols.  The  radiation  that  is  not  reflected 
or  scattered  and  reaches  the  surface  straight  forwardly  from  the  solar  disk  is  called  direct  or  beam  radiation.  The  scattered 
radiation  which  reaches  the  ground  is  called  diffuse  radiation.  Some  of  the  radiation  may  reach  a panel  after  reflection  from 
the  ground,  and  is  called  the  ground  reflected  irradiation.  The  total  radiation  consisting  of  these  three  components  is  called 
global  or  total  radiation.  There  are  three  components  of  the  solar  irradiance/solar  radiation  incident  on  inclined  plane  of 
slope  B and  orientation  [7] 

• The  direct  beam  component  described  by  the  Symbol  k. 

• The  sky  diffuse  component  described  by  the  symbol  Id. 

• The  ground  reflected  diffuse  component  described  by  the  symbol  Ir. 

HOURLY  BEAM  RADIATION  ON  INCLINED  SURFACE 

The  amount  of  beam  radiation  on  a tilted  surface  from  the  horizontal  surface  and  rotated  from  north  to  south  axis 
is  computed  by  multiplying  the  direct  horizontal  radiation  (Ib)  by  the  geometric  factor  (rb)  [8]. 


lb,  fi  =Ib  rb 

(1) 

cos  9 

IV  = 

D cos  Oz 

(2) 

Where  0z  is  the  zenith  angle  [22] 

cos0z  =sin  5 sin0  +cos  5 cos  0 cos  co  (3) 

5 is  the  declination,  0 the  geographic  latitude  and  co  is  the  hour  angle  calculated  at  midhour.  0 is  the  angle  of 
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incidence  for  an  arbitrarily  inclined  surface  oriented  toward  the  equator  and  calculated  by 

cos0  =sin  5 sin(0-P)  +cos  5 cos  (0-p)  cos  co  (4) 

HOURLY  GROUND-REFLECTED  RADIATION  ON  INCLINED  SURFACE 

Isotropic  Reflection 

Reflected  radiation  is  the  part  of  total  solar  radiation  that  is  reflected  by  the  surface  of  the  earth  and  by  any  other 
surface  intercepting  object  such  as  trees,  terrain  or  buildings  on  to  a surface  exposed  to  the  sky  is  termed  as  ground 
reflected  radiation.  A tilted  surface  at  slope  from  the  horizontal  has  a view  factor  Rr  to  the  ground  [18]. 

Rr  = (1-cosfi)  / 2 (5) 

Assuming  that  the  reflection  of  the  beam  and  diffuse  radiation  incident  on  the  ground  is  isotropic  and  that  the 
surroundings  have  a diffuse  reflectance  of  for  the  total  solar  radiation.  Subsequently,  the  reflected  radiation  (Ir)  from  the 
surroundings  on  the  surface  will  be  [21]: 

4=7/2  Ip  (1 -cos/3)  (6) 

where  Ir  is  the  ground-reflected  radiation  per  unit  area  arriving  at  the  inclined  plane,  4 is  the  angle  of  the  inclined 
surface,  I the  global  irradiance  on  a horizontal  surface  Wm2,  p is  the  ground  reflected  radiation  or  albedo  [9]. 

Anisotropic  Reflection 

Under  clean  and  cloudless  skies,  global  radiation  is  composed  primarily  of  direct  radiation.  When  the  ground  is 
covered  with  a layer  of  water  or  with  plants  having  glossy  leaves,  the  reflection  of  such  radiation  is  usually  anisotropic 
[10]. 

1=1/2  Ip (1 -cos P)[l+sin2  (02 /2)]{  \ cos  A| ) (7) 

Where  is  the  azimuth  of  the  tilted  surface  with  respect  to  that  of  the  sun? 

HOURLY  SKY  DIFFUSE  RADIATION  INCIDENT  ON  INCLINED  SURFACE 

Isotropic  Model  (Overcast  Sky  Condition) 

It  approximates  the  overcast  sky  condition;  the  sky  diffuse  radiation  incident  on  an  inclined  plane  is  given  by  the 
following  expression  [11]. 

Is  =1/2  Id(l  +cos  P)  (8) 

Where  Id  is  the  intensity  of  diffuse  radiation  onto  a horizontal  plane  W/m2. 

Anisotropic  Models  (Cloudless  Sky  Condition) 

The  anisotropic  model  is  applicable  only  to  a cloudless  sky  since  it  assumes  that  the  whole  diffuse  radiation 
comes  from  the  sky  sector  in  which  the  Sun’s  disc  is  located  at  a given  moment.  This  model  treats  diffuse  radiation  as 
beam  radiation  [12]. 

Ip,d  =Id^b  (9) 
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Anisotropic  Models  (Partly  Cloudy  Skies  Condition) 

Klucher  Model  this  model  for  conditions  of  partly  cloudy  skies.  This  model  is  based  on  a study  of  clear  sky 
conditions.  [8]  the  sky  diffuse  radiation  incident  on  inclined  plane  is  given  by  the  following  expression 

Is=l/2  Id  (1+  cos  p)  [1+F  sin3  (p/2)]  (1+F  cos2  0 sin3  0z)  (10) 

F=l-(Id  /I)2  (11) 

In  this  expression  the  parameter  F expresses  the  degree  of  anisotropy  as  a modulating  function  of  the  amount  of 
direct  radiation  received  by  the  surface. [13] 

TOTAL  HOURLY  RADIATION  ON  INCLINED  SURFACES 

The  total  amount  of  radiation  incident  on  an  inclined  plane  is  composed  of  beam,  ground-reflected,  and  sky- 
diffuse  components  [14].  At  locations  where  the  hourly  global  and  diffuse  radiation  on  horizontal  surfaces  are  known  or 
can  be  estimated,  the  global  radiation  on  an  inclined  surface  can  be  written 

I[iG  =Ip,b  + 2 Hr  +//i,  (12) 

This  equation  is  very  useful  in  analyzing  such  solar  energy  devices  as  flat  plate  Collectors.  [16-  18] 

METHODOLOGY 

Software  Tools  and  Model  Parameters 

The  hourly  solar  radiation  model  was  designed  by  using  MATLAB  program,  which  is  an  interactive  tool  for 
modeling  which  draw  a map  shows  the  temporal  and  spatial  distribution  of  the  solar  radiation  for  Iraq  and  the  neighboring 
countries.  Model  inputs  parameters  supplied  are  latitude,  longitude;  day  and  month  number,  solar  constant,  and  solar 
radiation  data  on  horizontal  surface  from  French  site  (SODA)  for  the  period  1/1/2005  to  31/12/2005.  And  in  this  work,  a 
constant  value  for  the  albedo  is  taken  equal  to  0.2  [23].  Figure  (1)  shows  the  flow  chart  of  simulation  program  for 
distribution  Hourly  solar  radiation  incident  on  inclined  surfaces  for  different  sky  condition. 
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Figure  1:  Research  Methodology 


RESULTS  AND  DISCUSSIONS 


The  hourly  global  solar  radiation  has  been  estimated  using  MATLAB  program,  which  estimates  hourly  global 
solar  radiation  using  isotropic  and  anisotropic  model.  The  method  presented  can  be  used  to  calculate  hourly  global  solar 
radiation  incident  on  inclined  surface  at  optimum  angle  has  been  calculated  at  different  sky  condition 
(Clear,  cloudy,  semi-cloudy). 

The  program  results  give  a clear  idea  about  the  behavior  of  hourly  solar  radiation.  It  is  revealed  from  the  results 
that  in  all  months  more  incident  solar  energy  radiation  on  tilted  surface  at  than  on  horizontal  surface  due  to  low  incidence 
angle  of  solar  radiation.  Following  Figures  present  the  simulation  results  obtained,  are  compares  Global  Solar  radiation 
Incident  on  horizontal  surface  and  Global  Solar  radiation  Incident  on  inclined  surface  at  different  sky  conditions.  In 
generally  the  monthly  variation  for  global  solar  radiation  for  Iraq  and  neighboring  countries  gives  us  that  the  highest 
intensity  value  is  in  the  summer  months  figure  (4),  while  the  lowest  is  in  the  winter  months  Figure  (2),  but  the  hourly 
variation  for  global  solar  radiation  during  the  morning  hours  (10:00-11:00)  we  notice  increasing  solar  radiation  in  all  the 
months  during  the  year,  In  other  words,  the  angle  inclination  to  be  effective  during  these  hours  it  was  observed  that  the  tilt 
angles  influence  maximum  hourly  energy  of  global  solar  radiation  especially  in  the  winter. 

Figures  (2,  6)  in  winter,  the  beam  solar  radiation  and  diffuse  incident  on  horizontal  surface  appear  close  to  each 
other  due  to  the  inclination  of  the  earth's  rotation  axis  at  northern  Hemisphere  away  from  the  sun,  which  makes  the  sun 
rays  incident  on  the  earth's  surface  at  very  small  acute  angles.  While,  in  summer,  the  beam  component  is  more  than  the 
diffuse  component  and  this  the  main  contribution  comes  from  the  beam  component. 
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Figure  2:  Temporal  and  Spatial  Distribution  for  Hourly  Global  Solar  Radiation  (W/m2)  Incident  on  Inclined 
Surface  for  Different  Sky  Condition  at  10:00-12:00  Local  Time  for  Day  15  in  winter  of  2005 
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Figure  3:  Temporal  and  Spatial  Distribution  for  Hourly  Global  Solar  Radiation  (W/m2)  Incident  on  Inclined 
Surface  for  Different  Sky  Condition  at  10:00  - 12:00  Local  Time  for  Day  15  in  spring  of  2005 
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Figure  4:  Temporal  and  Spatial  Distribution  for  Hourly  Global  Solar  Radiation  (W/m2)  Incident  on  Inclined 
Surface  for  Different  Sky  Condition  at  10:00  - 12:00  Local  Time  for  Day  15  in  Summer  of  2005 
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Figure  5:  Temporal  and  Spatial  Distribution  for  Hourly  Global  Solar  Radiation  (W/m2)  Incident  on  Inclined 
Surface  for  Different  Sky  Condition  at  10:00  - 12:00  Local  Time  for  Day  15  in  autumn  of  2005 


Figures  (6),  (7),  (8)  and  (9)  describe  the  distribution  of  solar  radiation  represent  on  different  width  lines  (25o-40o) 
considering  longitude  is  constant  (44o).  Although  the  increase  of  cloudiness  involves  the  reduction  of  direct  radiation  and 
the  rise  of  diffuse  radiation,  with  decrease  for  the  total  incident  energy,  but  when  comparing  between  these  figures  we 
notice  increasing  on  solar  radiation  during  the  morning  hours  and  after  the  mid  morning  while  the  surface  to  be  inclined  on 
optimum  angel  and  this  leads  to  increasing  on  number  of  hours  which  the  solar  radiation  become  high  intensive  at  different 
condition  sky  especially  at  Clear  sky  conditions  on  inclined  surface  Figure  (7). 

on  the  other  hand  the  hours  near  the  sunrise  and  sunsets,  the  Global  solar  radiation  values  are  low  because  the 
atmospheric  optics  mass,  zenith  angles  and  distance  are  increased  the  transferring  along  the  atmosphere  for  radiative 
fluxes,  which  make  the  occurrence  of  greater  attenuation  of  direct  radiation  possible  through  interaction  with  atmospheric 
constituents  by  the  processes  of  scattering,  absorption  and  reflection,  consequently  of  the  anisotropy,  and  in  the  afternoon, 
especially  in  winter  months,  average  hourly  energies  are  smaller  because  to  increase  of  cloudiness  and  the  accumulation  of 
water  vapor  during  the  day  Figure  (9). 


The  diurnal  evolution  is  dependent  of  photo  period  it  varies  10  am  and  13  pm  hours  for  the  winter  and  Autumn 
figure  (6)  (9)  and  varies  8 am  and  16  pm  hours  for  the  spring  and  summer  figure  (7)  (8)  because  application  of  angles  in 
different  sky  conditions  of  exposure  allows  the  lack  of  symmetry  in  the  diurnal  evolution  of  solar  radiation  when  compared 
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to  horizontal.  And  application  Liu  & Jordan  model  was  used  for  the  estimation  of  solar  radiation  incident  on  inclined 
surfaces  in  cloudy  skies  condition. 
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Figure  6:  Hourly  Global  Solar  Radiation  (W/m2)  at  Different  Latitude  in  winter 
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Figure  7:  Hourly  Global  Solar  Radiation  (W/m2)  at  Different  Latitude  in  spring 
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Figure  8:  Hourly  Global  Solar  Radiation  (W/m2)  at  Different  Latitude  in  summer 
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Figure  9:  Hourly  Global  Solar  Radiation  (W/m2)  at  Different  Latitude  in  autumn 
SPECIAL  CASE  FOR  BAGHDAD  CITY 


The  model  was  applied  to  estimate  the  global  radiation  on  inclined  surfaces  at  different  condition  sky,  using 
SODA  data  and  the  measured  data  incident  on  horizontal  surface  in  Atmosphere  Science  Center  in  Al-Mustansiriyah 
University  within  the  study  domain  by  Automatic  sun  tracker  with  pyrheliometer  at  longitude  44.2,  latitude  33.2  at  34 
meters  above  sea  level.  For  this  study,  one  year  for  solar  radiation  data  (Global  I and  Direct  lb)  for  year  2015.  Simulations 
were  compared  with  measurement  for  the  observational  data  for  solar  radiation  2015  are  available  Inclined  surface.  Figure 
(10)  Shown  similar  between  changes  of  solar  radiation  at  different  sky  condition  when  application  of  optimum  angles  for 
towards  the  equator.  And  shows  global  solar  radiation  at  local  noon.  It  was  observed  that  the  hourly  levels  of  global  solar 
radiation  incident  inclined  surface  (Ip5o)  are  superior  to  global  solar  radiation  incident  on  horizontal  surface  (IH,  g)  and  it 
indicates  that  the  difference  in  the  solar  energy  between  the  three  conditions  is  large. 
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(a)  (b) 

FigurelO:  (a)  Global  Solar  Radiation  (W/m2)  for  SODA  Data  in  9 Figure  (b) 
Global  Solar  Radiation  for  Measurement  Data 


CONCLUSIONS 

• It  was  found  from  this  study  the  values  increase  of  solar  energy  radiation  on  tilted  surfaces  in  all  months  for 
different  sky  condition  when  compared  with  the  values  of  solar  radiation  on  a horizontal  surface,  but  important 
case  in  cloudy  weather  the  surface  inclined  for  optimum  angle  and  this  cause  little  increasing  on  solar  radiation, 
This  is  important  in  the  months  (January,  February,  December). 

• In  the  morning  and  nightfall,  the  values  of  radiation  were  very  small,  but  between9  and  15  hours,  the  global  solar 
radiation  showed  the  highest  values  with  dependence  of  isotropy  and  anisotropy  radiation 

• When  applying  the  optimum  angle  of  inclination  cause  an  increase  in  the  sunshine  duration  during  the  winter 

• The  model's  ability  to  predict  and  describe  the  hourly  solar  radiation  during  different  times  and  different  condition 
sky. 

• Isotropic  and  anisotropic  models  predicted  more  incident  solar  energy  irradiation  on  tilted  surface  than  on 
horizontal  surface  due  to  low  incidence  angle  of  solar  radiation. 
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